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Abstract 

A new integral equation for three-magnon bound states in XXZ spin chain is 
suggested. Unlike the one presented by C. K. Majumdar about 40 years ago this 
equation does not have spurious solutions. Such an advantage is a result of de- 
composition of the wave function in the Bloch basis rather than in the basis of flat 
waves. 



1 Introduction 

An integral equation for three magnon bound states in Heisenberg ferromagnet was first 
suggested by C. K. Majumdar pQ and then directly obtained from spin Hamiltonian by 
several authors [2], [3]. The main advantage of this approach is its applicability to lattices 
with dimensions bigger than one. The main shortcoming is a presence of spurious solutions 
|3J,[1],[S] originated due to slightly inadequate representation for a three-magnon state 



|3) — M ni,n 2 ,n 3 S ni S n2 S n3 10). 
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Here |0) is the ferromagnetic ground state (S„ |0) = and the indices rij numerate a 
D-dimensional lattice. Amplitudes tt ni ,n 2 ,n 3 with n, = for j ^ k are unphysical and 
are separated in the Shrodinger equation from the physical ones [2], [3]. It is unphysical 
amplitudes that spurious solutions originate from. They have already been classified [6]. 
In the present paper for the case D = 1 we represent the three-magnon bound state in the 
form which contains only physical amplitudes and then obtain the corresponding integral 
equation. 



2 Hamiltonian and Shrodinger equation 

Hamiltonian of the S — 1/2 XX Z ferromagnetic chain 

H = - ^ ( S rt S n+l + S n S rtfl) + ^ ( S n S n+l - ~) . (2) 

n 

acts in an infinite tensor product of C 2 subspaces associated with each site of the chain. 
Here are the 5 = 1/2 spin operators corresponding to n-th site. For a three-magnon 
state we suggest the following representation 

|3>= £ e i (" 1+ " +p ) fc / 3 6 n _ m , p _„(fc)S m S,-S p |0). (3) 

m<n<p 

where 

io>=n®it>». w 

n 

is the ferromagnetic ground state. 

It is convenient to represent the corresponding Shrodinger equation in the form 

([H (k) + A(3 + V)]b(k)) = Eb m , n (k), (5) 

V / m,n 

where 

(H (k)b(kj) = ((S m ,i ~ l)6 m _i ( „(fe) + (5„,i - l)6 m+ i in _i(/c) - b mtn+1 (k)) 

\ / m,n Z V / 

Q—ik/3 , » 

(k) + (5 nj i - l)6 mjn _i(/c) - b 

m+l,n 

(Vb(k)) = -5 mi xb lj7l (k) - 5 n>1 b mil (k). (6) 
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(8) 



3 Three-magnon scattering states at A = 

In the free (A = 0) case Eq. (5) has the following system of /c-independent Bloch solutions 

^(0) /„\ _ \ ( Q j\qim+q 2 n\ _ Q i[qim+{q 1 ~q2)n} _|_ g[{.q2-qi)m-qw\ 
_ e *[(<?2-<?i)m+g 2 n] _|_ e -i[q2m+(q 2 -qi)n] _ g-ifem+gin] ^ q = ((^ q 2 ) (7) 

with dispersion 

E (k, q) = - cos Q - q^j - cos Q + q x - q^j - cos Q + g 2 
It is convenient to represent the wave function (7) in a compact form 

&£> n (q) = ^ V (-l)deg V[- (1) ^) m +- (2) fe)"] ; (_l)dcgc jfe! = (Q) 

where u give the following representations of £3 (the permutation group of three elements) 

^123(51, Q2) = (qi, 52), ^132(51, Q2) = (qi, 5i ~ 52), ^213(51, 92) = (52 - 5i, 52), 
W23i(5i, 52) = (52-5i,-5i), ^312(51,52) = (-52,5i -52), ^321(51,52) = (-52, -5i)- (10) 

From (9) follows the symmetry property 

^!>(<i)) = (-i) desw Oq), (n) 

which reduces f2, the fundamental region for the parameters gi and g 2 . If 

Q : < gi < 2tt, < g 2 < 2vr, (12) 

then 

Q = U ueSa u(n). (13) 

The system of functions (9) is normalized and complete 

00 

E ^(q)Oq) = %l-5~l)%2-5~2), (14) 
m,n=l 

3 = ^ r d qi r dqib^( q )b^M = ^s n ui5) 

{^) Jo Jo 
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In order to prove the normalization condition one should represent the product of wave 
functions in the sum (14) as follows 



1(0) (q)&(°) (g) = 1 (_l)dcg^+dcg [ i e -i[(^( 1 )( g i)-w( 1 )(gi))m+(^( 2 )( g2 )-<i( 2 )(g 2 ))r i ] 

= £(-l)«*»^ B (q-a;(q)), (16) 



where 



#mn(q) = - e -i[w (1) (gi)m+w( 2 )(g 2 )n] _ 1 ^g-i(<2im+<22n) _|_ e -i[qi(m+n)-q 2 n] 

_|_ g-i(g2m-gi(m+n)] _|_ e i[qim-q 2 (m+n)] _|_ ^[^(m+nj-gm] _|_ gifem+gin] ^ 



Since for each d m>n 

^ ^ (yd m ,m+n ~\~ <^ra+n,n) ^ ^ <^m,n ^ ^ ^n,n; (^) 
m,n>0 m,n>0 n>0 

one readily obtains from (17) 

1 oo oo 

9mA<l) = g[ E e^ im+ ^- E (e^ 1+ ^" + e^-^) + l 

m,n>0 m,n=— oo n=— oo 

= ^[47r 2 % 1 )% 2 )-27r(% 1 + g 2 ) + 5(g 1 -g 2 )) + l] . (19) 
Hence according to (16) and (19) 

E C(q)e*(q) = ^ E (-i) dcsH %i - ^ (1) (<zi))% 2 - *< 2 >(&)), (20) 

m,n>0 w€53 

(other terms cancelled after averaging over 5 3 ). But w(q) G f2 only for w = wi 23 . So all 
terms in (20) with u 7^ Wi 23 should be omitted and we get Eq. (14). 

In order to prove completeness let us first define an action of S 3 on the configuration 
space by the following formula 

w (1) (gi)m + u (2) (q 2 )n = giw (1) (m) + q 2 uj {2) {n). (21) 

According to (10) and (21) 

wi 23 (m, n) = (to, n), wi 32 (to, n) — (to - ra, — n), w 2 i 3 (m, rc) = (— to + n), 
w 23 i(to, n) — (-to - n, to), w 3 i 2 (to, n) = (n, -to - n), w 32 i(to, n) = (— n, -to). (22) 
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In this notations 

&- 0) -(q)& (0) (q) = - ^2 (-l) degw+deg V [9l(w(1)(m) ~^ (1)( ™ ))+92(w(2)( ™ ) "' : ' (2)( ™ ))] (23) 



6 



So 

r2n /•2tt 



7535/ ^S!n(q)^!n(q)=E(- 1 ) deS ^^ 1 )M^)(n)- (24) 

Since m,n,m,n > Eq. (15) follows from (22) and (24). 

4 Shrodinger equation in the q-space 

Let 

b m ,n(k) = ^J n b( nU^ b (k, q)dq, (25) 

or equivalently 

6(fc,q)= ^ 5{2n(q)6™,»(*). ( 26 ) 

m,n>0 

be a decomposition of the wave function in the basis of Bloch functions. In this framework 
Eq. (5) takes the form 

(E (k, q) + 3A - E)b(k, q) = -6A / V(q, q)6(fc, q)dq, (27) 

where 

CO 

^(q,q) = -^E^^M^+SSC^q))- (28) 

n=l 

As it follows from (11), (25) and (28) 

6(fc,w(q)) = (-l)** w &(*,q), VMq),q) = ^(q, W (q)) = (q, q). (29) 

Hence the integration over Q, in Eq. (27) may be extended to integration over tt Q 

(E (k, q) + 3A - E(k))b{k, q) = -A / dq, / dq 2 V(q, q)6(fc, q), q G ft . (30) 

Jo Jo 

For scattering states 

q) = £ (M) + 3A, (31) 
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and hence the right side of (30) should be zero. For bound states we assume 

k 

E bound (k) < min^Escauik, q)) = 3A - 3 cos -. (32) 
Let us now calculate V(q, q). According to (16) and (28) 

1 oo 

V(<L q) = E (- 1 ) dC ^ E 9U^) - q) + 5n,i(^(q) - q)- (33) 

we53 n=l 

As it follows from (17) 

E^(q) + fn,i(q) = g E [ e ~ % ( e ~* 92n + e i<?2(n+1) + e i(M)n + e i(91 - g2)(n+1) ) 

n=l n=l 

_|_ e »(?l-?2) ^g*?!™ _|_ e -%(™+ 1 ) _|_ e -iQ2n _|_ e 92(n+l) j _|_ ^2 ^ e «(<?2-<?i)n _|_ e i(qi-?2)(n+l) 

+e mn + e -^(n+i)^j = I ^ e^ (1) ^)5(a; (2) (g2)) - ^(cosgi + cosg 2 + cos ( 9l - ft)) .(34) 

we 53 

Eqs. (33) and (34) result in 

^(q,q) = --i- V" (-l) de s-- de s-e l [- (1) ^)-- (1) ^)](5[a;( 2 )(g 2 ) -cl; (2) (g 2 )]. (35) 

127T ^— ' 

A substitution of (35) into (30) and evaluation of the sum over u using appropriate 
changes of q gives for a bound state 

(E (k,q) + 3A-E bound (k))b(k,q) = ^ T dq ± T dq 2 Y^(-l) dc ^ 

.e i [- (1) ^)-5i] 5 [ a; (2)( ?2 ) _ g ~ 2 ] 6(A;?q ) ; (36 ) 

or equivalently 

(E (k, q) + 3A - E bound (k))b(k, q) = A E (-l) dcg ^ (1) ^V(A;, ^ 2 )(g 2 )), (37) 

we<S3 

where 2 

^(fc,5) = J- [*dqer*b(k,q,q). (38) 

Eqs. (37) and (38) result in an integral equation 

«m - £ E (-')«" f ft(t , , J" w *) (39) 

tJfc03 
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or in an extended form 

2vr J E (k,q,q) + 3A - E bound (k) LV / 
+ ^(9-2,-) _ g-itf+aj^^ _^ + ^ e -*(9+fl _ l^k, q - g)] . (40) 

For future simplication we notice that 

g, g) + 3A - £ 6ound (*0 = 7(M(* ~ *+)(* ~ *-) ? ^ 

z \ z - ~ z +) 

where 

~ _ iq _ 3A - E bound (k) - cos (fc/3 + g)± jjje, g) , . 

2 - e ' z ± - (e<(*/3-9) + e- lfc / 3 ) ' 1 j 

and 



7 (fc,g) = y(cos(- + g)+ J E; f)oW (A;)-3A) - |e«*/3-g) + e -^| . (43) 
From (42), (43) and (32) follows that 

z + z^ = z + z_ = 1, \z+\ > 1, \z-\ < 1. (44) 



Hence from (41) and (44) one readily gets 

"2* (!_ e %-2g))^ l-e^ 



or 



(45) 



2vr 7 E Q (k,q,q) + 3A - E bound (k) j(k,q) 
Now using (45) and the following formulas 

#o(fc, -5, ?) = E (k, q + q,q) = E (k, -q, q), (46) 

one may reduce (40) to the form 

/ A(l - e iq z 2 _ )\ A r 2n [e i( -i+ 2 V - + e -%~+2<?) _ 

V 7(M )^ ,9 '~2^Jo E (k,-q, q) + 3A - E bound (k) 



(47) 



1 _ A(l-e^) \ A /^ [cos 3(q + g)/2 - cos (g - g)/%(fc, gdg 

7(*>9) / t 7o #o(&, -9> ?) + 3A - E bound (k) 

where 

<p(k,q) = e i *' 2 i/>(k,q). (49) 

Eq. (48) is the main result of the paper. Although it is analogous to Eq. (91) of 
Ref. 1, its derivation does not need introduction of any additional constructions (as it 
has been for the Majumdar equation [2], [3]). In the Bloch basis Eq. (48) directly follows 
from representation (3). According to (26), (38) and (49) up to a normalization constant 

oo 

<p(k, q)=Y, M^ ( " +1/2) + &i, n (£)e^ ( " +1/2) . (50) 
n=l 
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5 Checking on the Bethe Ansatz result 

An exact form of the three- magnon wave function is well known [TJ. Namely 

b min (k) = z?- l (k)z%- l (k), (51) 



where 



zi(k) = ^-^Ae^ + e" 2 ^ 3 ), z 2 (k) = Zl (k) = -^i_(2Ae- lfc / 3 + e 2jfc / 3 ). (52) 
In fact it may be readily proved that (51) satisfies the system (6) with 



2 V z 2 (k) z 1 (k)J 2 V z 2 (k) 

8A 3 + cosfc 

3A ~ 4A»-1 ■ (53) 



According to (50) and (51) 



e «//2 e -ig/2 



9) = : 77T + 77T- ( 54 ) 

Since 

-g, g~) + 3A - £ 6oMnd (fc) = 7(fc ' g) _ ( f_" ~ ^ , (55) 

z^z. — z + j 

Eq. (48) takes the form 

Uk,q) + A^z 2 - 1)) ( r— — -r-r + . e " W/ ' 1 = A( ^; ~ * +) I (e 3 ^ 2 z 2 
\' K ' v ') \e- l i - zx(k) e*« - z 2 {k) J 2m J\ 

-e-^~z - e-/ 2 + e" 3 -/ 2 i) ( f + 1 ) _ ^f. (56) 

z/ \1 — zzi(/cj z(z — Z2(/c)J / (z — z + )(z — z_) 

Integral in the right side of (56) may be decomposed and readily calculated 

Z l - ' / ,3i 9 /2~3 _ -« ? /2~2 _ iq/2~ , -3iq/2\ <lz 



2m' J V / (1 — zz\{k)){z — z + )(z — z_) 

= (e^/2^3 _ e - lq /2-2_ _ Qlq /2-_ + e - 3lq /2^ 1^ (57) 

v 1 — Z-Z\[k) 

and (w = 1/z) 

!_1± I ( e »9/2 5 _ e -^ 2 - e^ 2 - + e^ 2 ^ -—^ r , 

2m J \ z z l J {z - Z2{k)){Z — z + )(z — Z-) 

= E^Zl± I ( e 3^/2 _ ^ q /2 w _ e i q /2 2 + e -M q /2 w 3} dw 

2ni J v ' (1 - wz 2 (£;))(l - z+iw)(l - z_w) 

3 3V2 _ e -n/2]_ _ e i q /2 1 + e -%/ 2 l N j !+ (5g) 

2+ 4 Z+/ 2+ - Z 2 (fc) 



Using (57), (58) and (44) one may reduce Eq. (56) to the form 




_ e- 3 ^/ 2 ^^ - l)( e 2iq z- - 1) 

1 — Z-Zi(k) 

e 3iq/2 (e~ i(1 z 2 _ - l)(e- 2iq z_ - 1) 
1 - z-z 2 (k) 



,(59) 



which may be checked directly by rather cumbersome calculation. 

6 Conclusions 

In the present paper we obtained an integral equation (48) for three-magnon bound states 
in ID Heisenberg ferromagnet. The suggested equation is based on the representation 
(3) which does not contain unphysical amplitudes related to spurious solutions. The 
derivation is based on the decomposition (25) of the wave function in the basis of Bloch 
wave functions. The obtained equation was directly tested on the Bethe Ansatz solution 
(51). Basing on this result we suggest that for a study of bound states in a complex model 
it is better to decompose a wave function not in the flat waves basis but in a basis of 
states related to an appropriate solvable model. 

The author is very grateful to P. P. Kulish for a formulation of the problem and to M. 
I. Vyazovsky for an interest to the work. 
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